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We study multipartite entanglement in non-equilibrium quantum phase transition (NEQPT) attainable in a 
coherently driven atomic ensemble undergoing collective decay. The NEQPT arises in the steady state of the 
system as the drive field strength is varied in comparison with the decay parameter. A close connection is ob- 
served between the critical behavior and the multi-partite entanglement in the system determined via the von 
Neumann entropy. The derivative of the von Neumann entropy shows A -type behavior typical of the second or- 
der phase transitions. We further show that the calculated bi- and tri- partite entropies satisfy the Lieb inequality. 
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Quantum phase transitions (QPT) — phase transitions oc- 
curring in quantum matter at zero temperature — signify criti- 
cal behavior in the ground state of the system under the in- 
fluence of an external physical parameter flJL£| and are dif- 
ficult to study via direct methods due to large number of 
degrees of freedom involved. Recently, there has been a 
strong interest in the simulation of such complex systems 
(e.g., condensed-matter systems) in simple controllable quan- 
tum systems that are simpler to understand, calculate or ob- 
serve such as trapped-ions [2]. Some notable examples of sim- 
ulation of familiar QPTs of condensed-matter systems include 
the observation of the superfluid to Mott insulator transition in 
Bose-Einstein Condensates (H and study of the BCS to BEC 
crossover in degenerate Fermi systems [4]. In this context, 
study of QPTs, for example as observed in spin chains [5], 
which could be simulated in Dicke type systems [6] or optical 
lattices |7[], is especially important to understand the connec- 
tion between entanglement and the critical behavior. It has 
also been noted that the universality of the critical phenom- 
ena simply extends to the pairwise-entanglement Hunl- Fur- 
ther, there have been extensive studies of the so-called Lip- 
kin model involving chain of fermions with short-range inter- 
actions lldl ; however, optical realizations of this model are 
far from clear. Study of the spin chains is further impor- 
tant because of possible applications in quantum communi- 
cation [11]. It could, however, be noted that so far there has 
been no study relating the multipartite entanglement and crit- 
ical behavior in non-equilibrium QPT (NEQPT). 

Here we show how a non-equilibrium quantum phase tran- 
sition (NEQPT) can be achieved in a collective atomic en- 
semble and study the relation of multiparticle entanglement 
with this type of critical behavior. The dynamics of the atomic 
system can be controlled by the external drive field, variation 
of which over a wide-range, in comparison with the collec- 
tive decay parameter, gives rise to a phase transition in the 
steady state of the system. The system consists of a collec- 
tively driven ensemble of two (or three) level atoms that is 
simultaneously subjected to a collective emission of radiation 
as for example induced by a cavity. Thus, due to the con- 



tinuous driving and collective decay the system is inherently 
a non-equilibrium system. The non-equilibrium steady state 
(NESS) density matrix of the system can be determined and 
allows study of the entanglement of the complete system or 
parts of it in terms of the von Neumann entropy which has 
been recognized as a measure of entanglement and quantum 
correlations for optical systems for quite some time |T2, 131. 

We consider a collection of two level atoms, described by 
the spin-half operators (Sj,S^) and the transition frequency CO, 
interacting resonantly with a single mode of the cavity field, 
described by boson operators a and af. The atoms are in addi- 
tion driven by a coherent laser light resonant with the atomic 
transition. The Rabi frequency of the interaction of the atom 
with the coherent light is |i2|e"^ and with the cavity mode is 
g. The total Hamiltonian is given by 

=ha>Y,S Z j + ftffl^a - I'Y.^^ S j a + H - c 

j 

-h\Q.\ 1 £(^e- ia>, Sfe i 1'j + H. c.) (1) 

The dynamics described by the above Hamiltonian is to be 
supplemented by the dissipative terms describing the leakage 
from mirrors. The cavity field decay rate is taken to be k. In 
the bad cavity limit, K ^> g, cavity field can be adiabatically 
eliminated to arrive at the following master equation for the 
atomic density matrix: 

p = i \Q.\[e- i 4>S + +e i ' t 'S-,p} -F{S + S-p -2S + pS~ + pS + S~) . 

(2) 

Here F — \g\ 2 / K is the collective decay parameter, and 5 ,± = 
^jSje 1 ^ are the collective atomic operators. Here the den- 
sity matrix is written in a frame rotating with the frequency 
of the applied field CO therefore it does not appear explicitly 
in the dynamical equation. We add that another very popu- 
lar system, three-level system with single-photon emission on 
the \a) — \b) transition, which has been found to be useful for 
long-distance quantum communication 1141, could also be de- 
scribed by the same dynamical equation <(2j (See Fig.[T](B)). 

The collective atomic operators are analogous to the ladder 
operators for a spin N /2 system consisting of total N two-level 



atoms contained in the ensemble. It is important to note at 
this point that there are two competing dynamical processes 
going on in the system arising due to the external coherent 
drive and the collective decay. As it will be shown later this 
competition gives rise to the phase transition in the system 
as the relative amplitudes of |fi| and T are varied. We have 
recently showed how to generate generalized Werner states, 
which are analogues of the two-particle maximally-entangled 
mixed states, by coherently driving the collectively decaying 
ensemble of two-level atoms IU5I1 . 

Despite the non-equilibrium nature of the system under 
consideration there exists a steady state that can be formally 
written as llrjll 

p = (R-)-\R + )- l /D=(S-+iG)- l (S + -iG)- l /D, (3) 

where R = S~ + i(\Q.\/T)e^, D = Tt[(S~ +iG)- 1 (S+ - 
iG) ] is the normalization and the newly introduced param- 
eter G is given by G = iQ./F. This allows us to determine 
the steady state of the system, nevertheless, it is important to 
know the typical time required to attain this steady state in or- 
der to determine the rate at which the system can be swept 
through the phase transition. In general, the time required 
for such a system of N atoms to arrive at the steady state is 
roughly Vt « l/N for any arbitrary value of Q. (See Ref. 11511 
for discussion of a two atom case). Thus, the time-rate of the 
variation of the drive field Rabi frequency, required to scan 
through the phase transition, needs to be smaller than the time- 
rate of approach to the steady state. Within this regime Q. in 
Eq. d2J is taken to be time-independent. The collective de- 
cay rate increases as the size of the system increases allowing 
faster variation of the amplitude \£l\ to study the passage of 
the system through the phase transition. 

In the spin representation of the collective system, the 
steady state density matrix can be written as: 

P = T> t HG)- m (iG*)- n (S-) m (S + y (4) 
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Representing the spin operators in Eq. dU in the basis of 
angular momentum states with S = N/2 (so called Dicke 
states [17]) the individual density matrix elements p p>q in 
Eq. (|5]l can be determined. This is possible as the dynamics 
in Eq. (f2]i is invariant under the exchange of the particle labels 
and the state remains completely symmetric and is confined 
to the manifold with S = N/2 provided the initial state of the 
atoms is completely symmetric. These symmetry considera- 
tions are also important in the extraction of density matrix of 
the different-sized parts of the system for further study of the 
entanglement. 

Once the density matrix for the complete system is 
known, the density matrix p^ M > for part of the system 
containing M particles can be projected out via p' M ) = 



(A) 

|a>" 



|b>. 



(B) 

|a>- 



FIG. 1: (A) Level scheme for two-level atoms: A coherent drive 
and collective decay into the cavity mode allow realization of the 
dynamics in Eq. (f2}. (B) Level scheme for three-level atoms: A co- 
herent drive and microwave field on the lower transition alongwith 
collective decay into quantum vacuum of a cavity allow for a dif- 
ferent realization of the dynamics in Eq. {2]l that has been also use- 
ful for quantum communication 1 14], In the notation of this figure 
st = \a){b\. 



The projected density matrix can then be used to study 
properties of a part of the system. To maintain the symmetry 
of the parts of the system, the state of the N particle state has 
to be expressed as a superposition of completely symmetric 
parts with M and N — M particles. This is possible through 
the spin representation that we have already introduced; thus, 
\N/2,p) = E p=;+ft C^f |M/2,j)|(iV-M)/2,fc) , where the 

C N -^ are the Clebsch-Gordan coefficients arising due to the 
coupling of the component angular momenta to form the total 
spin N /2 system 111 811 . Thus, the density matrix elements 



pj$ = ( f , M | P (M) | f , V) are given by 

Nil <N-M)/2 

p»= I p m E E E<Xf- 

p,q=-N/2 l=-(N-M)/2 p=H+l q=v+l 

(6) 

allowing determination of multipartite entanglement in the 
system. In particular, we study bi- and tri-partite entangle- 
ment in an atomic ensemble containing total N particles. 

In Figs. |20 we show a connection between bipartite en- 
tanglement and the critical behavior in the NEQPT observ- 
able in the collectively driven system as a function of the 
scaled Rabi frequency associated with the drive field Cl s = 
2G/N = 2i £1/ 'NT . The plot in Fig. [2] shows the von Neumann 
Entropy of the entangled state determined from the projected 
two-particle density matrix p' 2 ' as [19] 



S(p( 2 )) = -Tr(p( 2 )log 2 p( 2 )). 



(7) 



-M/2 
~l±.y=-M/2 
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It is known that the non-equilibrium phase transition occurs in 
these systems at £l s — 2G/N = 1 Illa,l20fl ; we observe that the 
entanglement character of the two-partite entanglement in the 
Af particle system follows the phase transition. In the asymp- 
totic limit of N, G — > °° with D. s finite we observe a sharp phase 
transition in the entropy and correspondingly the bipartite en- 
tanglement in the system at £l s = 1 . The phase transition char- 
acteristic can be seen in an accentuated manner via the plot of 
derivative of the von Neumann Entropy of the projected bi- 
partite system in Fig. [3] From the plot it can be noted that the 
nature of variation of the derivative of the Entropy through a 
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FIG. 2: Bipartite entanglement in the non-equilibrium QPT exam- 
ined via the von Neumann Entropy defined in Eq. 0. 



FIG. 5: Tripartite entanglement in the non-equilibrium QPT exam- 
ined via the von Neumann Entropy S(p' 3 )) = — Trp( 3 )log 2 




FIG. 3: Bipartite entanglement in the non-equilibrium QPT exam- 
ined via the derivative of the von Neumann Entopy. The behavior in 
the asymptotic limit is similar to the A -type variation of the specific 
heat typical of second order phase transitions. 



phase transition is similar to the A -type variation of the spe- 
cific heat observed in second order phase transitions. 

We next discuss the role of the phase of the driving field. In 
fact, we demonstrate that coherent nature of the drive is sig- 
nificant in leading to phase transition behavior. In this case 
the non-equilibrium steady state density matrix averaged over 
the random phase possesses only the diagonal elements that 
are non-zero. The averaged von Neumann entropy of the pro- 
jected system is plotted in Fig. [4] It can be observed that the 
transition, even in the asymptotic limit, is not as sharp as in the 
case of coherent drive and the entropy increases more-or-less 
in a smooth manner. 

From the study of bipartite entanglement it can also be 
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FIG. 4: Bipartite entanglement in the non-equilibrium QPT exam- 
ined via the average of the von Neumann Entropy over random 
phases of the driving field. There is no sharp transition even in the 
asymptotic limit. 



observed that an average two-particle entanglement in the 
asymptotic limit for large £l s is very close in value to the bi- 
partite entanglement in the two-particle system giving a clear 
indication that the two particle entanglement is not a suffi- 
cient measure of entanglement in the N particle system. It 
is furthermore important to note that the popular measure of 
bipartite entanglement, so-called, concurrence II2U \22 J. is not 
a suitable measure for the study of phase transitions |g]; our 
conclusions are similar and we therefore do not discuss it here. 

It is, of course, important to study the multipartite entangle- 
ment in the system as it undergoes NEQPT. The techniques we 
have presented here can be employed to study entanglement 
of any part, containing M particles, of the total A^-partite sys- 
tem. We present the results of this study below in the form 
of the tripartite entanglement. We note in this case, however, 
that it is difficult to determine the the asymptotic behavior. 

In Fig. [5] we study tripartite entanglement via the von Neu- 
mann Entropy S(p( 3 >) — — Trp( 3 'log 2 p' 3 ) contained in the 
NESS of the N particle system. The trend of entropy with 
the increasing value of the coupling parameter £l s = 2G/N is 
the same as observed for bipartite entanglement in Fig. |2]ex- 
cept for the sharper features near the transition point £l s = 1 
and larger value of the overall tripartite entropy. Thus, the 
result is similar as in the context of spin chains where the en- 
tropy of the multiparticle part increases with size in a logarith- 
mic fashion ]8[] . It can be noted that the asymptotic value of 
the von Neumann entropy of the M partite system in the limit 
£l s — > oo is given by log 2 (M+ 1) where the state of the system 
approaches to a completely mixed state. The asymptotic value 
of the entropy for bi- and tri- partite systems turns out to be 
1 .585 and 2 respectively and is reconfirmed by the plots. 

It is, furthermore, important to study true multiparticle en- 
tanglement by separating out entanglement in smaller parts of 
it. This can be accomplished for example via something called 
relative entropy. The relative entropy of a projected M particle 
density matrix p ( M ) can be evaluated as 



M 

s*(p (m) ) = E 

k=l 



M 



(M-k) 



(8) 



where S(p( M " fc )) is the von Neumann Entropy obtained from 
Eq. for the density matrix of M — k particles projected out 
of the total N particles. The Relative entropy of tripartite en- 
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FIG. 6: Tripartite entanglement in the non-equilibrium QPT exam- 
ined via the von Neumann Relative Entropy defined as 3S(p' 2 )) — 
3S(pO) — S(p( 3 )) where S is defined in the caption of Fig. [5] 
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FIG. 7: Diagnostic of the determined von Neumann Entropy of 
the parts of the system via the Lieb inequility S(pO) +S(pw) < 
2S(p( 2 '). The negative values of the quantity plotted show that the 
Lieb inequality is satisfied. 



tanglement is shown in Fig. [6] The trend is similar to the 
complete tripartite entanglement, that is with increasing N the 
transition happens closer to the critical point £2 S = 1 signi- 
fying close connection between entanglement and the critical 
behavior of the NEQPT. 

As a diagnostic measure of the calculated results we check 
that the von Neumann entropy of the projected parts of the 
system satisfies the Lieb inequality [19] in Fig. [7] The result 
is as expected and the quantity S(p^) +5(p^) -2S(p (2) ) 
remains negative for the whole parameter range. 

To summarize, we show realization of a non-equilibrium 
quantum phase transition and its connection with multipar- 
tite entanglement in a collectively driven atomic ensemble 
that is simultaneously subjected to collective decay. This 
study has a lot to offer: First of all it allows one to study 
non-equilibrium quantum phase transitions in simple control- 
lable systems and their associated entanglement properties. It 
also offers possibility to study correspondence between meso- 
scopic and macroscopic systems undergoing QPTs. In view 
of our results it would also be worthwhile to study the first 
order nonequilibrium phase transitions as in optical bistable 



systems 02311 from the point of view of quantum entanglement. 

GSA thanks the National Science Foundation, grant no 
nsf-ccf 0524673, for financial support. 



* Electronic address: KT-Kapale@wiu.edu 
^ Electronic address: Girish.Agarwal@okstate.edu 
[1] S. Sachdev, Quantum Phase Transitions (Cambridge University 

Press, Cambridge, UK, 1999); P. Coleman and A. J. Schofield, 

Nature 433, 226 (2005). 
[2] D. J. Wineland et al, Physica Scripta T76, 147 (1998); D. 

Leibfried et al, Phys. Rev. Lett. 89, 247901 (2002); T. Schaetz 

et al, Appl. Phys. B 79, 979 (2004). 
[3] M. Greiner et al, Nature (London) 415, 39 (2002). 
[4] J. T. Stewart, J. P. Gaebler, C. A. Regal, and D. S. Jin, Phys. 

Rev. Lett. 97, 220406 (2006). 
[5] X. Wang and P. Zanardi, Phys. Lett. A 301, 1 (2002). 
[6] S. Schneider and G. J. Milburn, Phys. Rev. A 65, 042107 

(2002); N. Lambert, C. Emary, and T. Brandes, Phys. Rev. Lett. 

92, 073602 (2004); C. F. Lee and N. F. Johnson, Phys. Rev. 

Lett. 93, 083001 (2004); J. Reslen, L. Quiroga, and N. F. John- 
son, Europhys. Lett. 69, 8 (2005); V. Buzek, M. Orszag and M. 

Rosko, Phys. Rev. Lett. 94, 163601 (2005). 
[7] J. J. Garcia-Ripoll, M. A. Martin-Delgado, and J. I. Cirac, Phys. 

Rev. Lett. 93, 250405 (2004). 
[8] A. Osterloh, L. Amico, G. Falci, and R. Fazip, Nature (London) 

416, 608 (2002). 
[9] V. E. Korepin, Phys. Rev. Lett. 92, 096402 (2004). 
[10] A. Garg, Europhys. Lett. 22, 205 (1993); S. Dusuel and J. Vi- 

dal, Phys. Rev. Lett. 93, 237204 (2004); R. G. Unanyan and 

M. Fleischhauer, Phys. Rev. Lett. 95, 260604 (2005); R. G. 

Unanyan, C. Ionescu, and M. Fleischhauer, Phys. Rev. A 72, 

022326 (2005). 

[1 1] J. S. Pratt and J. .H. Eberly, Phys. Rev. B 64, 195314 (2001); S. 

Bose, Phys. Rev. Lett. 91, 207901 (2003). 
[12] H. Huang and G. S. Agarwal, Phys. Rev. A 49, 52 (1994). 
[13] S. M. Barnett and S. J. D. Phoenix, Phys. Rev. A 40, 2404 

(1989). 

[14] L.-M. Duan, M. D. Lukin, J. I. Cirac and P. Zoller, Nature (Lon- 
don) 414, 413 (2001); L.-M. Duan and H. J. Kimble, Phys. Rev. 
Lett. 92, 127902 (2004). 

[15] G. S. Agarwal and K. T. Kapale, Phys. Rev. A 73, 022315 
(2006). 

[16] R. R. Puri and S. V. Lawande, Physica 101A, 599 (1980). 
[17] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[18] I. I. Sobelman, Atomic Spectra and Radiative Transitions, 2nd 

ed. ( Springer- Verlag, Berlin, 1996). 
[19] M. A. Nielsen and I. L. Chuang, Quantum computation and 

quantum information (Cambridge University Press, Cambridge, 

2000). 

[20] L. M. Narducci, D. H. Feng, R. Gilmore, and G. S. Agarwal, 
Phys. Rev. A 18, 1571 (1978). 

[21] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 

[22] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[23] L. Lugiato, in Progress in Optics, edited by E. Wolf (North- 
Holland, Amsterdam, 1984), Vol. XXI, pp. 69-216. 



